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Electronic structures of organic molecule encapsulated BN nanotubes under 

transverse electric field 
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The electronic structures of boron nitride nanotubes (BNNTs) doped by different organic 
molecules under a transverse electric field were investigated via first-principles calculations. The 
external field reduces the energy gap of BNNT, thus makes the molecular bands closer to the BNNT 
band edges and enhances the charge transfers between BNNT and molecules. The effects of the 
electric field direction on the band structure are negligible. The electric field shielding effect of 
BNNT to the inside organic molecules is discussed. Organic molecule doping strongly modifies the 
optical property of BNNT, and the absorption edge is red-shifted under static transverse electric 
field. 

PACS numbers: 75.75.-l-a, 73.22.-f, 72.80.Rj 
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I. INTRODUCTION 

Nanotube is a very important actor in nanodevice ap- 
plications due to its novel properties. Compared to the 
more widely studied carbon nanotube (CNT), boron ni- 
tride nanotube (BNNT) has its own distinctive proper- 
ties. Besides the excellent mechanical stiffness and ther- 
mal conductivity comparable to CNT^i BNNT has a uni- 
fied electronic structure regardless of its diameter, chiral- 
ity, and the number of walls of the tube*?''^-^ BNNT also 
shows pronounced resistance to oxidation, and it is sta- 
ble up to 700 °C in air, while CNT survives only below 
400 °C.^ AU these properties make BNNT an attractive 
candidate for nano-electronics. 

Pristine BNNT is a wide-gap semiconductor. For elec- 
tronics application, it is desirable to make p- or n-type 
doping to BNNT. Our previous study shows that the 
electronic structure of BNNT can be modified by or- 
ganic molecule encapsulation i^ Electrophilic molecule in- 
troduces acceptor states in the wide gap of BNNT close 
to the valence band edge (VBE), which makes the doped 
system a p-type semiconductor. However, with typical 
nucleophilic organic molecules, instead of shallow donor 
states, only deep occupied molecular states are observed. 
There is a significant electron transfer from BNNT to an 
electrophilic molecule, while the charge transfer between 
a nucleophilic molecule and BNNT is negligible. 

On the other hand, previous theoretical studiesii^i^ 
showed that the band gap of CNT or BNNT can be re- 
duced and even closed by applying a transverse electric 
field, due to the giant Stark effect (GSE).^ The electric 
field can mix the nearby subband states in the valence 
band (VB) complex and conduction band (CB) complex 
separately, leading to an electric field dependence for the 
band gap. This effect is more remarkable in BNNTs than 
in CNTs, due to the reduced screening of the electric 
field in BNNTs. The GSE was lately confirmed exper- 
imentally using the bias dependent scanning tunneling 
microscopy and scanning tunneling spectroscopy J^ 

Since the electronic structure of BNNT can be tuned 



by both organic molecule encapsulation and transverse 
electric field, it is thus interesting to see what will happen 
if we apply both. In this paper, by performing density 
functional theory (DPT) calculations, we study the elec- 
tronic structures of organic molecule encapsulated BN- 
NTs under a transverse electric field. The electrostatic 
shielding of BNNT to the inside molecules and the ef- 
fects of organic molecule doping and electric field on the 
optical properties of BNNT are also studied. 

PoUowing our previous study^S several typical elec- 
trophilic and nucleophilic molecules are considered 
for BNNT encapsulation. The two electrophilic 

molecules studied are tetracyano-p-quinodimethane 
(TCNQ) and tetrafluorocyano-p-quinodimethane 
(F4TQ). Three nucleophilic molecules are selected: 
tetrakis(dimethylamino)ethylcnc (TDAE), anthracene 
(ANTR), and tetrathia-fulvalene (TTF). The (16,0) 
BNNT is chosen as a prototype in this study. We name 
the organic molecules capsulated BNNT as M@BN, 
where M is the name of the doped organic molecules. 




TCNQ@BN 
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FIG. 1: (Color online) Optimized geometrical structures of 
the MQBN systems under a 0.5 V/A electric field. B in dark- 
salmon, N in blue, C in grey, H in yellow, F in aqua, and S in 
goldenrod. 



The rest of the paper is organized as follows: In Sec. II, 
the theoretical approach and computational details are 
briefly described. In Sec. Ill, the calculated band struc- 
tures, charge transfers between BNNT and molecules, 
shielding effect under electric field, and optical proper- 
ties are discussed. Finally, we conclude in Sec. IV. 



II. COMPUTATIONAL DETAILS 

To investigate the geometry structure and electronic 
states of the organic molecule encapsulated BNNT, we 
performed first-principles spin-polarized DFT calcula- 
tions. The computational details of geometry optimiza- 
tion and band structure calculation can be found in 
our previous study,-£ and we only give a brief summary 
here. We used the projector augmented-wave (PAW) 
metho d^^i^^ implemented in the Vienna Ab Initio Sim- 
ulation Package (VASP)J^ii^ The Perdew, Burke and 
Ernzcrhof (PEE) exchange-correlation functional was 
chosen. ^^ We note that more rigorous band structure 
calculations are required to yield correct band gapsi-^^ 
However, PBE functional has been widely used in similar 
systems, and it should give correct trends which are in- 
terested in this study. The energy cutoff for plane- waves 
was 400 eV. In our calculation, each BNNT was separated 
by 10 A of vacuum, and the minimum distance between 
periodically repeated organic molecules along the BNNT 
tube axis was ~8.7A. A 4x1x1 F-centered Monhkorst- 
Pack fc-point grid was used for Brillouin zone sampling. 

The way to handle electric field in VASP is adding an 
artificial dipole sheet in the middle of the vacuum part 
in the periodic supercell.^'^'^* In our study, the direction 
of the electric field was chosen to be perpendicular to the 
nanotube's axis (along the y direction), with a magni- 
tude of 0.20 or 0.50 V/A. As discussed later, the selected 
fields are strong enough to make significant change of the 
electronic structure, and weak enough to avoid artificial 
field emission. 



Kramer-Kronig transformation, and the imaginary part 
£2 can be calculated based on the independent-particle 
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FIG. 2: Electrostatic potential (in volt) across the BNNT in 
the M@BN systems under 0.5 V/A external electric field. The 
top of the highest occupied bands are marked by dash lines. 
The numbers on the top of the curves are the effective electric 
field (in V/A), calculated by the slope of the line connecting 
the left high-potential point and the right low-potential point 
in vacuum. 



All linear optical properties, such as absorption co- 
efficient and refractive index, can be obtained as func- 
tions of dielectric function. The real part of the dielec- 
tric function can be obtained from its imaginary part by 



approximationjii i.e.. 

All 



UJ 



(1) 



where f2 is the unit cell volume, and lu is the photon 
energy. Summations for i and j are performed for va- 
lence and conduction bands, respectively. P^ denotes 
the dipole transition matrix elements obtained from the 
self-consistent band structure calculations and Wu is the 
k point weighting. The same method has been applied to 
h-BN, and the results are in reasonably good agreement 
with experiments.^" The number of bands used for di- 
electric function calculations are the sum of the electron 
number and the ion number. The number of occupied 
bands is half of the number of electron, and the num- 
ber of unoccupied bands are more than 380 for all the 
systems studied in this paper. The optical polarization 
in dielectric function calculations is perpendicular to the 
nanotubc axis. In this way, for some parts of the tube 
wall, the optical polarization is nearly perpendicular to 
the BN layer, while for the other parts of the tube wall, 
it is roughly parallel to the BN layer. 



III. RESULTS AND DISCUSSION 
A. Geometric and band structures 

As reported in our previous study;^ after the organic 
molecules are encapsulated into BNNT, their geometries 
changes little. This suggests the weak interaction be- 
tween BNNT and the encapsulated molecules, due to the 
large distances between them. The separations between 
neighboring molecules and between molecule and BNNT 
are larger than additions of the corresponding atomic 
van der Waals radius. Optimizations under the trans- 
verse electric field gave similar structure as in the zero- 
field case. The optimized geometrical structures under 
a 0.5 V/A electric field are shown in Fig. [TJ The re- 
laxed geometries of the M@BN systems before and after 
applying electric field change little. The largest change 
of bond length, comparing to the geometry under zero 
field, is smaller than 0.01 A. Consistently, the direction 
of the transverse field does not affect the geometric and 
electronic structures. Our test calculation gave almost 
identical band structures for F4TQ@BN under 0.5 V/A 
transverse electric field along two different directions. For 
simplicity, in the rest part of the paper, all electronic 
structure calculations were done based on the geometri- 
cal structures obtained without electric field. 

Before calculating the electronic structure of M@BN 
with transverse electric field applied, two important is- 
sues are needed to be discussed. One is that the vacuum 
region should be broad enough to avoid the overlap be- 
tween the charge density and the artificial dipole sheet 



used to generate electric field. The width of the vacuum 
layer between two neighboring BNNTs in our calcula- 
tion is 10 A, and the work function of (16,0) BNNT is 
6.36 eV. The charge density of BNNT roughly decay as 
exp(— v^?"), where cf) is the work function in eV and r is 
the distance from the BNNT in Ai^i The charge density 
thus decays one magnitude every 0.91 A. The dipole sheet 
is 5 A far from the BNNT, and there is thus only neg- 
ligible charge density there. The other issue is that the 
vacuum region should not be too broad so as to attract 
electrons to vacuum (field emission). In Fig. [51 we plot 
the electrostatic potential of the M@BN system under a 
0.5 V/A electric field, the strongest electric field used in 
our calculations. We can see that the lowest potential in 
vacuum region is still higher than the edge of the high- 
est occupied band (HOBand), which confirms that the 
artificial field emission has not been introduced in this 
study. 
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FIG. 3: Band structures (from F to X) of the pristine BNNT 
under (a) 0.0, (b) 0.2, and (c) 0.5 V/A electric fields. Fermi 
levels are marked by dashed lines. Density distribution for 
HOBand, LUBand, and LUBand-f-l in pristine BNNT under 
(d-f) 0.0 and (g-i) 0.5 V/A electric fields. The electric field 
direction is from up to down. 

We first calculate the band structure of the pristine 



BNNT under the 0.2 and 0.5 V/A electric fields, and the 
results are shown in Fig. [31 The band gap decreases from 
4.45 to 3.73 and 1.98 eV, when the field increases from 
to 0.2 and 0.5 V/A, respectively. Our result is consistent 
with the previous theoretical study on similar systems 4 
Such band gap decrease is a result of the extensive mix- 
ing among the subband states within the valence and 
conduction band complex under electric filed. The band 
mixing is also reflected by the band-resolved charge den- 
sity distributions, as shown in Fig. [31 where the density 
of the HOBand, the lowest unoccupied band (LUBand), 
and LUBand+1 under 0.0 and 0.5 V/A field are plotted. 
The orbital density move along the electric field for occu- 
pied states, and against the field direction for unoccupied 
states. From the density distributions, the LUBand state 
can be identified as a nearly free electron (NFE) stated 
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FIG. 4: Band structures (from F to X) of the pristine BNNT 
and the MQBN systems under 0.5 V/A electric field. Fermi 
levels are marked by dashed lines. 

The band structures of the M@BN systems under 0.2 
and 0.5 V/A electric field are then calculated. Similar 
to results under zero electric field,— the bands are ei- 
ther mainly from molecules or mainly from BNNT. The 
molecular bands are much flatter than the BNNT bands. 
Since there is no qualitative difference between the band 
structure of the M@BN systems under 0.2 and 0.5 V/A 
electric fields, we will focus on the 0.5 V/A case (Fig. [J). 
For the two electrophihc molecules (TCNQ and F4TQ), 
the lowest unoccupied molecular bands become very close 
to the BNNT VBE, which means good p-type doping. It 
also suggests a considerable charge transfer between the 
molecules and the BNNT. We define the gap between the 
lowest unoccupied molecular band and the BNNT VBE 
as Ep, and the gap between the highest occupied molec- 
ular band and the BNNT conduction band edge (CBE) 
as En- Their values are listed in Table [H For the nu- 
cleophilic organic molecules TDAE and TTF, En drops 
from 2.30 and 2.85 eV at zero field^. to 1.08 and 1.54 
eV, respectively. However, it is still too large to form 
good n-type doping. For ANTR@BN, the highest occu- 
pied molecular band is pushed into the BNNT valence 
band manifest, and there is only an unoccupied molecu- 
lar band close to the BNNT CBE. Therefore, neither n- 
type or p-type doping is introduced to BNNT by ANTR 



encapsulation. 

To consider the molecular effect on the intrinsic elec- 
tronic structure of BNNT under electric field, we define 
the energy gap between the CBE and VBE of BNNT 
as E^^ . As shown in Table HI without external field, 
it is not very sensitive to the type of molecules. How- 
ever, E^'^ for TCNQQBN and F4TQ@BN under elec- 
tric field is much larger than other systems under the 
same field. This can be easily understood by their elec- 
tronic structures. Under strong external field, the energy 
gaps of TCNQ@BN and F4TQ@BN become very small, 
which means enhanced metallicity and stronger electro- 
static screening. With stronger screening, the reduction 
of electrostatic potential along the direction of electric 
field is slower (see Fig. [2|), therefore, E^^ reduction 
compared to pristine BNNT and other M@BN systems 
is smaller. Similar behavior has been observed in our 
study on defective BNNTs under an electric fieldi^ 



B. Charge transfer 



LUBand has a high density (i?^ < 3.5 A)i^ 
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FIG. 5: (Color online) The cylinder-integrated differential 
electron density curves for the the pristine BNNT and the 
M@BN systems under a 0.5 V/A electric field. The blue verti- 
cal lines indicate the radial boundary between molecules and 
the BNNT, in the viewpoint of charge transfer. The black 
vertical line in the first panel marks the radius of the BNNT. 



An important property of the molecule encapsulated 
system is the charge transfer between the molecules 
and BNNT. We calculate the charge transfer using the 
method described in our previous workf^ originally pro- 
posed by Lu et al.ji^ Briefly speaking, the boundary be- 
tween the organic molecule and the BNNT is determined 
by the maximum/minimum position of the cylinder- 
integrated differential electron density curve, and the 
corresponding maximum/minimum value is the charge 
transfer value. The differential electron density is defined 
as the difference between electron densities of the M@BN 
system under electric field and the sum of charge densities 
of BNNT and molecules without electric field. The range 
of the radial position i?b of boundary between molecule 
and BNNT can be estimated by the density distribu- 
tion of the BNNT bands. For electrophilic molecules, Rb 
should be inside the region where the BNNT HOBand 
has a high density {Rjj < 6 A), while for nucleophilic 
molecules, Rb should be inside region where the BNNT 



For pristine BNNT, the differential electron density is 
just the density change before and after applying elec- 
tric field. As shown in Fig. [SJ the electron density in- 
creases inside the nanotube wall, while it decrease outside 
the wall. Of course, the absolute value of this change is 
very small (~0.01 electron). For the electrophilic organic 
molecules TCNQ and F4TQ, under 0.5 V/A electric field, 
there are 0.28 and 0.38 electrons per molecule transfer 
from BNNT, about twice as the value without field (0.10 
and 0.24 electron).^ This is due to the closer molecular 
state in energy to VBE. The boundary radius is 4.3 A, 
which is basically the midpoint between the molecule and 
the BNNT as in the no field case. For the nucleophilic 
organic molecules TDAE, TTF and ANTR, the charge 
transfers form molecules to BN also increase, which is 
0.006, 0.011, and 0.008 electron per molecule, respec- 
tively (0.004, 0.008, and 0.005 electrons without external 
field) . Such small electron transfers also indicate that the 
doped system is not a good n-type semiconductor. 



TABLE I: Energy gaps Ep, E„, and Eg^ of the pristine 
BNNT and the M@BN systems under electric field £. The 
energy gaps are in eV, and the electric field £ is in V/A. 



Gap 


£ BN 


TCNQ F4TQ 


TDAE TTF ANTR 


Ep 


0.0 — 


0.15 


0.09 


4.37 


3.55 


2.96 




0.2 — 


0.17 


0.11 


— 


3.25 


2.65 




0.5 — 


0.09 


0.05 


— 


— 


1.93 


En 


0.0 — 


— 


— 


2.30 


2.85 


3.71 




0.2 — 


— 


— 


1.86 


2.49 


3.40 




0.5 — 


— 


— 


1.08 


1.54 


— 


E^ 


0.0 4.45 


4.45 


4.48 


4.45 


4.46 


4.43 




0.2 3.73 


4.00 


4.18 


3.63 


3.79 


3.81 




0.5 1.98 


2.79 


2.91 


2.31 


2.14 


2.18 



C. Shielding effect 

It is interesting to study the effects of BNNT to the 
geometry and electronic structure of the inside molecules 
under an electric field. In Fig. [SI we compare the molec- 
ular partial densities of states (PDOS) with and with- 
out external field for two representative molecules: an 
electrophilic one (F4TQ) and a nucleophilic one (TTF). 
The calculated PDOS shows different behaviors for these 
two types of molecules. For the electrophilic organic 
molecules F4TQ (Fig. [5^), there is an nearly rigid up- 
shift of the molecular PDOS after applying the electric 
field. This is due to the electron transfer from BNNT 
to the molecules. For the nucleophilic organic molecules 
TTF (Fig. [Bt), the molecular PDOS almost does not 



change after applying the electric field, which is consis- 
tent with the little charge transfers between BNNT and 
the molecules. 
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FIG. 6: (Color online) Molecular PDOS of (a) F4TQ@BN and 
(c) TTFQBN under 0.0 (solid) and 0.5 (dashed) V/A electric 
field. DOS of individual (b) F4TQ and (d) TTF molecules 
under 0.0 (solid) and 0.5 (dashed) V/A electric field. 



We also check the shielding effect by applying electric 
field to individual molecules without a BNNT encapsu- 
lation. First, we relax the geometry of the individual 
molecules under the 0.5 V/A electric field. Contrasting 
to the encapsulated case, we found a significant geometry 
relaxation after applying the electric field. Then, we fix 
the molecular geometry to the one as inside BNNT, and 
calculate the electronic structure under zero and finite 
electric field. The effect of electric field to the molecular 
electronic structures for individual molecules is different 
from that for molecules encapsulated by BNNT (see Fig. 
in]). For F4TQ molecule, the DOS energy shift caused by 
electric field is smaller than that for PDOS of F4TQ@BN. 
However, we can see a notable change of the overall DOS 
shape after applying the electric field. For TTF, elec- 
tric field shifts the DOS of individual molecules, but it 
does not change the molecular PDOS. In summary, our 
results indicates, although not good as CNT, the BNNT 



still provides a significant electrostatic shielding for the 
inside molecules: the geometry of the inside molecules 
and the electronic structures of nucleophilic molecules 
are almost not affected by an external field under the 
protection of BNNT. 



D. Optical properties 

The calculated imaginary part of the dielectric func- 
tion of the pristine BNNT without external electric field 
is shown in the inset of Fig. [7^. A much denser fc-point 
sampling is required to obtain smooth curves. ^S However, 
the main features of €2 are already presented in our re- 
sults. And our calculated €2 of prinstine BNNT agrees 
well with the previous theoretical result on the zigzag 
BNNTs.?° Mainly, there are three peaks at about 5.6, 
10.6, and 13.5 eV, respectively. The strong peak around 
5.6 eV corresponds to the tt — *■ tt* interband transitions. 
The two other peaks higher in energy and weaker in in- 
tensity are associated with interband transitions also in- 
volving the a bands. We note that the peak around 10.6 
eV does not appear in our test calculation with the opti- 
cal polarization along the tube axis. By optical selection 
rule analysis, ^° the 10.6 eV peak should correspond to 
TT ^ cr* and cr — > TT* transitions, and the 13.5 peak mainly 
contributed by cr —> a* transitions. In Fig. [7^, we also 
plot £2 of the pristine BNNT under transverse electric 
field. Similar to a previous study on pristine BNNT,— 
the whole imaginary dielectric function is not strongly 
affected by the static field, but the absorption edge is 
red-shifted. 
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FIG. 7: (Color online) The imaginary part (52) of the theoret- 
ical dielectric function of the BNNT and the M@BN systems, 
with (red) and without (blue) the external static electric field. 
Inset: £2 plotted in a larger energy scale. 



With molecule intercalation, all the three peaks are 
still clearly shown in the M@BN system. However, new 
features appear in the low-energy region from to 5 eV. 



These peaks below the BNNT absorption edge (~4.5 eV) 
are clearly from transitions between molecular bands and 
the BNNT bands or between molecular bands themselves. 
For the two electrophilic molecules, there are peaks close 
to the zero energy. This is consistent with the very nar- 
row band gap calculated for these two systems. There 
is no significant peaks below 4.5 eV for TDAE@BN, and 
the new peaks for TTF@BN and ANTR@BN are close 
to the BNNT absorption edge. Our results indicate that 
the optical property of BNNT can be strongly modified 
by different organic molecule encapsulation. 

Next, we study the case where both molecule interca- 
lation and electric field exist. The theoretical imaginary 
part of the dielectric functions of the molecule doped 
BNNT under a 0.5 V/A transverse static electric field are 
also presented in Fig. [T] The three peaks from BNNT 
are still not strongly affected by the static field. How- 
ever, if we look into the details below 5 eV, we can see 
three different behaviors. For transitions between valence 
band and conductance bands of BNNT, the electric field 
GSE broadens the band manifests, thus broadens the ab- 
sorption peaks too. This broadening effect red-shifts the 
absorption edge of BNNT as shown in Fig. [7Ji. For 
transitions between molecular bands, due to electrostatic 
shielding of BNNT, the electric field does not strongly af- 
fect the peaks, as shown in Fig. ^ and[7f. For transitions 
between BNNT bands and molecular bands, the electric 
field affects both the amplitude of the peaks and the peak 
positions, as shown in Fig. [Tls and[71:. 



IV. CONCLUSION 

In summary, we calculate the electronic structure and 
optical properties of BNNT under transverse electric field 



with different molecule intercalations. When applying 
transverse electric field, the band gap of BNNT decrease, 
while the band structures of organic molecules change 
little, or only with a rigid shift. After applying electric 
field, there is more charge transfer between molecules and 
nanotube in the p-type semiconductors TCNQ@BN and 
F4TQ@BN. BNNT provides relatively good electrostatic 
shielding for the geometry relaxation of the inside organic 
molecular chain, and for the electronic structure of nu- 
cleophilic molecules. Organic molecule doping strongly 
modifies the optical properties of the composite M@BN 
systems. The absorption edge is slightly red-shifted un- 
der static electric field. 
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